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Abstract 


Recently, several ultraluminous X-ray (ULX) sources were shown to host a neutron star (NS) accretor. We perform 
a suite of evolutionary calculations, which show that, in fact, NSs are the dominant type of ULX accretor. 
Although black holes (BH) dominate early epochs after the star-formation burst, NSs outweigh them after a few 
100 Myr and may appear as late as a few gigayears after the end of the star-formation episode. If star formation is a 
prolonged and continuous event (i.e., not a relatively short burst), NS accretors dominate the ULX population at 
any time in the solar metallicity environment, whereas BH accretors dominate when the metallicity is sub-solar. 
Our results show a very clear (and testable) relation between the companion/donor evolutionary stage and the age 
of the system. A typical NSULX consists of a~ 1.3 Mz NS and ~1.0 Mz Red Giant. A typical BH ULX consists of 
an8 Mo BH and ~6 M., main-sequence star. Additionally, we find that the very luminous ULXs (Ly = 104! 
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ergs 


') are predominantly BH systems (~9 M3) with Hertzsprung- gap, donors (~2 M. 


). Nevertheless, some 


NSULX systems may also reach extremely high X-ray luminosities (>10*! erg s~'). 
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1. Introduction 


Ultraluminous X-ray source (ULX; for a review, see Feng & 
Soria 2011) is defined by two observational properties: 


1. it is a point-like (i.e., not extended) off-nuclear X-ray 
source with a peak emission localized in the X-ray band; 

2. it emits isotropic equivalent X-ray luminosity in excess of 
10°° erg s~!, which is approximately the Eddington limit 
(EL) for a spherically accreting stellar-mass black hole 
(sMBH; ~10 Mo). 


Until recently the main problem with unveiling the nature of 
the ULXs stemmed from the absence of reliable measurements 
of the masses of the accreting compact objects. This has 
changed with the discovery of pulsing ultraluminous X-ray 
sources (PULXs), which contain neutron stars (NS). The first 
X-ray pulsar was discovered in M82 X-2 (Bachetti et al. 2014). 
Then, two other PULXs have been identified: P13 in NGC 
7793 (First et al. 2016; Israel et al. 2017b) and NGC5907 
ULX1, (First et al. 2017; Israel et al. 2017a). All of these 
PULXs show regular pulses with periods ~1 s, characteristic of 
NS accretors. 

These discoveries proved that some ULXs do contain stellar- 
mass compact objects. Moreover, it is possible that NSs may 
reside in a vast majority of the ULXs (NSULXs; Kluzniak & 
Lasota 2015; King & Lasota 2016; King et al. 2017). Existence 
of PULXs provides strong evidence in support of the models 
involving a super-Eddington accretion onto a compact star 
(e.g., Begelman 2002; Poutanen et al. 2007; King 2009) as an 
explanation of the ULXs phenomenon (e.g., Kording et al. 
2002). 

However, apparent super-Eddington luminosities can also be 
reached without breaching the EL if the radiation of an 
accreting compact star is beamed (e.g., King et al. 2001; 
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Poutanen et al. 2007). Detailed accretion disk simulations 
appear to support the importance of beaming (e.g., Ohsuga 
et al. 2005; Ohsuga & Mineshige 2011; Jiang et al. 2014; 
Sadowski & Narayan 2015). 

In this study, we perform massive simulations in order to 
uncover the nature of the stars forming the ULX population. 
Motivated by the recent observational progress, we limited our 
investigation to the stellar-mass compact objects with super- 
Eddington and/or beamed emission. 


2. Simulations 


The calculations were performed with the use of the 
StarTrack population synthesis code (Belczynski et al. 
2002, 2008) with significant updates described in Dominik 
et al. (2012) and Wiktorowicz et al. (2014). 

We start the evolution of each system from the zero-age 
main sequence (ZAMS), which we assume to happen at the 
same moment for both stars. The primary has an initial mass of 
M, = 6-150 Me drawn from a power-law distribution with 
index —2.7 (Kroupa & Weidner 2003). The mass of the 
secondary (M,) is chosen from the 0.08 to 150 Mz, range to 
preserve a uniform mass-ratio distribution. Typical predeces- 
sors of X-ray binaries (XRB) with NS or BH accretors should 
reside in such a range of masses. The initial separations are 
uniform in logarithm (P(a) ~ 1/a; Abt 1983). The eccentri- 
cities’ distribution is the thermal-equilibrium (Duquennoy & 
Mayor 1991). The binary fraction was set to 50% for stars with 
ZAMS mass below 10 Mj and 100% for systems with more 
massive primaries. We are aware of the results of Sana et al. 
(2012) who provide other relations for initial parameter 
distributions, but their results were obtained for stars with 
masses 15-60 Mz and were limited to solar metallicity. 
Therefore, we decided to keep the previously established 
distributions. For comparison between our adopted and Sana 
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Figure 1. Schematic representation of the dominant NSULX evolutionary 
routes in young (top; Rns,ms) and old (middle; Rns,rq) stellar populations. 
Bottom: typical BHULX evolution (SF regions). The age and mass values 
reflect those in a typical system. For an explanation of abbreviations, see the 
text (Section 3.1). For phases in which the mass of one of the components 
changes rapidly, two values are provided: the mass at the beginning of a phase 
and at the end (in parentheses). 


et al. (2012) initial distributions the reader may refer to de 
Mink & Belczynski (2015). 

We simulated the evolution of 2 x 10’ binary systems for 
every model and scaled the results to the Milky Way equivalent 
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Figure 2. Time evolution of the number of ULXs since the beginning of 
star formation. The BHULXs (thick line) appear early, but 100-1000 Myr 
later the ULX population becomes dominated by the NSULXs (dashed 
line), except constant SF models with sub-solar metallicities. Number of 
ULXs for a star-forming mass of 6 x 10'° Mz (see Section 2). Presented 
values are for all ULXs (including less typical routes, which are not 
presented in Table 1). For the relative abundance of young and old ULXs, 
see Figure 5. 
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Figure 3. Similar to Figure 1, but for the evolution of the most luminous 
(Lx > 10*! erg s~!) ULXs. Top/bottom plot is for a BH/NS accretor. 


galaxy (Mmwec = 6 x 10!°M, (Licquia & Newman 2015); 
star-formation rate (SFR) equal 6.0 Mz yr~! for constant star 
formation (SF) during 10 Gyr, and 600 M; yr7! for burst SF 
with a duration of 100 Myr). We do our study for solar 
metallicity (Zo), 10% of solar (Z>/10), and 1% of 
solar (Z,/100). 


2.1. Accretion Model 


We implemented a model based on nonlinear scaling of X-ray 
luminosity and mass accretion rate (eg, Shakura & 
Sunyaev 1973), 


Lx tot = Leaa(1 + In mr), (1) 


where 7 is the mass transfer (MT) rate in Eddington rate units. 
Mass accretion onto the compact object (Myce) is limited to 
Meaa and the rest of the mass transferred to the system is lost in 
a wind from the inner disk region. Such outflows from ULXs 
were recently observed by Pinto et al. (2016). 


2.2. Beaming Model 


The collimation of radiation from the innermost parts of the 
accretion disk, i.e., beaming, has a prodigious influence on the 
apparent isotropic luminosity. The beaming factor is defined as 
b = Q/4n, where 2. is the solid angle of emission (e.g., King 
et al. 2001). If we assume that there are two conical beams with 
opening angle 6, then Q = 47(1 — cos@/2). 
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The isotropic equivalent X-ray luminosity (Lx) can be 
expressed as 


Lx tot 
Lx i (2) 
Under the assumption of the isotropic distribution of disk 
inclinations in space (see King 2009), the probability of 
observing a source along the beam is equal b. 

The results of theoretical studies as well as the outcomes of 
detailed modeling of accretion disks favor the dependence of 
beaming on the MT rate (e.g., Lasota 2016). King (2009) 
proposed that b scales as 


bol pets 
m 


b=1 m<85. (3) 


For very high MT rates, this prescription may give an extremely 
small value of b, i.e., strong beaming. For example, a Galactic 
ULX candidate $S433 has ri ~ 3000-104 (e.g., Fabrika 2004) 
for which Equations (3) and (2) give b 7 x 10-7-8 x 10-6 
and Lx = 10*6 erg s-!. Therefore, we will assume saturation at 
m = 150. For higher mass accretion rates, we will consider that 
the beaming is constant and equal b = 3.2 x 10-3(@ = 9°). 
Theoretical background for beaming saturation was presented by 
Lasota et al. (2016). 


3. Results 


In this section we provide results for the most typical systems 
and focus on the standard model (Section 2). In the Appendix, 
we provide additional characteristics of the simulated population 
of ULXs and results for other models. 


3.1. Formation of NSULXs 


According to our results, NSULX are present in ULX 
populations of all ages and metallicities. The first NSULXs form 
as early as ~6 Myr after the start of SF, but the ULX phase may 
occur also in a very old system (fage ~ 5 Gyr). Dominant 
evolutionary routes leading to the occurrence of high MT rate 
and formation of an NSULX depend on the age of the system at 
the time of the ULX phase. Below, we describe the typical routes 
for early (Rns.ms), mid-age (Rng HG), and old (Rys.RG) stellar 
populations (see also Figure |). Ranges represent 50% of values in 
present day populations.’ 

6-59 Myr; Route Rnsms: Typical companion MS—On 
ZAMS, the masses of the primary and secondary are 
7.2-12 Mz and 1.3—2.0 Mo, respectively. The initial separation, 
500-2600 Ro, shrinks during the CE phase commenced by the 
primary. After the CE phase, the primary forms an NS through 
a core collapse supernova explosion (SN). The natal kick 
leaves the system on an orbit that is tight enough for the 
secondary to fill its Roche lobe (RL) within about 2 kyr due to 
nuclear expansion and gravitational radiation (GR) before it 
leaves the main sequence (MS). During the RLOF, the 
companion is comparable in mass (0.9-1.5 M3) to the NS 
(.1-1.4 Mz), thus the MT is stable and may be sustained as 
long as ~1.3 Myr (see Figure 1, top). 

400-820 Myr; Route Rys.uc: Typical companion HG—In 
approximately 54%-73% percent of the systems following the 


7 For example, populations as we predict them to look for an observer able to 
see all the objects at the same time. Problems related to light propagation (e.g., 
light speed, interstellar absorption) are not discussed in this study. 
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Table 1 
Typical Parameters of ULXs 
Present ZAMS 
Route 1(Myr) A t(Myr) M(M>2) M,(Mo) a(Ro) Mzams,a(Mo) Mzams,b(Mo) azams(Ro) 
Z=Z- 

Reu.MsS 11-22 <0.2 7.7-8.6 5.9-7.2 18-22 40-50 6.6-12 3800-4600 

ev. route: CE1(4-1; 7/8-1) SN1 MT2(14-1) 
Rns,MS 6-38 S0.4 1.1-1.3 1.1-1.4 2.2-3.6 8.2-9.3 1.3-1.7 700-1200 

ev. route: CE1(3/4/5-1; 7/8-1) SN1 MT2(13-1) 
RNs,HG 430-730 0.3-0.8 ~1.3 0.6-1.0 7.3-17 7.9-8.3 1.6-2.0 800-1200 

ev. route: CE1(3/4/5-1; 7/8-1) SN1 MT2(13-1/2) 
FRNS,RG 2300-4300 0.1-0.2 ~1.3 ~1.0 35-50 7.5-1.9 1.2-1.5 1100-1700 

ev. route: CE1(6-1; 12-1) MT2(12-3) AICNS1 MT2(13-3) 

Z = Z_/10 
RBu.MS 4-37 <27 5.7-8.9 5.6-7.8 14-19 26-37 6.0-11 1800-3700 
Rns.Ms 9-55 S0.8 1.1-1.3 1.2-1.5 1.6-2.7 7.2-12 1.5-1.9 1000-2600 
RNs.HG 400-820 S0.6 1.2-1.3 0.6-0.9 9.3-22 6-8 1.3-1.8 600-1500 
FRNS,RG 1400-2600 S0.1 ~13 ~1.0 35-45 6.3-6.8 1.3-1.6 1700-2700 
Z = Z5/100 

RBu.MS 8.7-19 <1.0 12-21 7-11 12-17 35-55 9.1-14 540-2000 
Rs,Ms 15-59 0.6-1.3 1.2-1.4 0.9-1.0 2.8-3.5 7.2-8.2 1.7-2.0 500-900 
RNs.HG 400-690 S0.4 ~13 0.6-0.8 17-35 6.0-6.3 1.7-2.1 700-900 
RNS,RG 1200-2100 S0.1 ~13 ~1.0 30-45 6.0-6.2 1.3-1.7 900-1300 


Note. Typical present and ZAMS parameters of most common ULX evolutionary routes. Ranges represent 50% of values in present day populations. The parameter 
columns represent, respectively: age of the system in ULX phase (), duration of ULX phase (Ar), compact object mass (M,), donor mass (M,), separation (a), primary 
mass on ZAMS (Mzams,a), secondary mass on ZAMS (Mzams,b), Separation on ZAMS (azams). The schematic evolutionary routes provided for Z = Zs are the same 
for other metallicities. Symbols represent (Belczynski et al. 2008) CEl—common envelope (donor: primary); MT2—mass transfer (donor: secondary); SN1— 


supernova (primary); AICNS 1—accretion induced collapse (primary); 1 


MS; 2—HG; 3—RG; 4—CHeB; 5 


evHeS; 12—ONeWD; 13—NS; 14—BH. 


early AGB; 6—thermal-pulsing AGB; 7—HeS; 8— 


Table 2 
Typical Parameters of HLXs 
Present ZAMS 
Route t(Myr) A t(Myr) M,(M5) MyM 5) a(Ro) Mzams,a(Mo) Mzams.p(Mo) azams(Ro) 
Z=Zo 
ReBuUG 15-40 S0.08 8.0-10 1.7-3.3 20-90 40-50 3.8-8.2 3700-4600 
ev. route CE1(4/5-1; 7/8-1) SN1 MT2(14-1/2) 
RNs.evHes 17-40 <S0.001 1.3-1.4 1.7-2.6 <44 10-11 8.5-10 30-300 
ev. route MT1(2/3/4/5/6-1/2/4) SN1 CE2(13-3/4/5; 13-7/8) MT2(13-8/9) 
Z=Z,5/10 
RBuH.HG 20-60 <S0.06 7.2-9.5 1.7-3.3 12-120 25-35 4.0-7.0 2000-5000 
RNS, evHes 40-60 <0.001 ~13 1.5-2.0 8.4-13 8.8-9.5 5.5-7.3 170-430 
Z = Z,/100 
RBuG 20-45 <S0.04 10-18 1.2-3.7 10-75 25-45 5.0-9.0 700-2600 
FNS. evHes 50-75 <0.001 ~13 1.9-2.5 2.1-6.9 6.4-6.8 5.5-6.2 170-280 


Note. Typical present and ZAMS parameters of typical most luminous ULXs (HLX; routes Rgy yo and 7Rns,evHes). The table is organized in the same way as Table 1. 


Rns.s route, the RLOF during the MS phase, if present, is not 
strong enough to power a ULX. These secondaries start to 
reach their terminal-age MS and expand rapidly several 
hundred Myr after the start of the SF. If the separation is 
short enough, they will fill their RL. Thus, at this time, a 
Hertzsprung-gap (HG) star becomes the most common NSULX 
companion. 

1200-4300 Myr; Route Rus rg: Typical companion RG— 
The evolution of the late-time NSULXs is significantly 
different than that of the earlier ULXs. After a CE phase, the 
primary forms an oxygen—neon white dwarf (ONeWD) instead 


of an NS, and thus the mass of the primary is on average lower 
compared with routes Rys.ms and Rys Ho. However, when the 
secondary becomes a red giant (RG) and fills its RL of 
~ 10-20 Ro, the primary accretes additional mass and forms an 
NS due to an accretion induced collapse (AIC). Afterward, the 
RG refills its RL and a short ULX phase occurs (Figure 1, 


middle). 


3.1.1. Formation of BHULXs 


In general, BHULXs are a minority in our results. Never- 
theless, they dominate the ULX population during the initial 
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Table 3 

Frequently Used Abbreviations 
ULX Ultraluminous X-ray source 
BHULX ULX with a BH accretor 
NSULX ULX with an NS accretor 
PULX Pulsing ULX 
MS Main sequence 
HG H-rich Hertzsprung gap 
RG Red giant 
CHeB Core helium burning 
HeS Helium star 
evHeS Evolved helium star 
SN Supernova 
AIC Accretion induced collapse 
RL Roche lobe 
RLOF Roche lobe overFlow 
MT Mass transfer 
GR Gravitational radiation 
SF Star formation 
WD White dwarf 
ONeWD Oxygene-Neon WD 


phase of the constant SF (the exception being populations with 
low (<Z,/10) metallicity), and during the SF burst. They are 
also more numerous than NSULXs among the ULXs with 
luminosity in excess of 10*' erg s~' (see Sections 3.2 and 3.3). 

4-37 Myr; Route Rgyms: Typical companion MS—The 
initial binary on ZAMS consist of a massive primary 
(26-55 M3) and an intermediate-mass secondary (6-14 M,). 
The initial semimajor axis is ~540-4000 Ro. 

The heavier star evolves very quickly and during the CHeB 
phase fills its RL, commencing the CE episode. As a result, the 
primary is deprived of its hydrogen envelope and the separation 
drops to a few x 10 Ro. A direct collapse occurs shortly 
afterward. The combined action of the companion’s nuclear 
expansion and GR leads to an RLOF. The phase of ULX 
emission is short, S1.0 Myr (Figure 1, bottom). 

In the case of the gums channel, it is essential that the 
heavy primary expands and fills an RL of several x100 Ro, 
which terminates the expansion. There may be some observa- 
tional evidence that for the most massive stars (VM > 40 M3) 
the large dimensions are avoided (Mennekens & Vanbeve- 
ren 2014). Additionally, the massive stars are almost all found 
in close interacting binaries (Sana et al. 2012), where very 
quickly their envelope is removed by RLOF and/or CE. It is 
naturally expected in our models. Since the post-MS expansion 
is very rapid (~200 kyr), and envelope removal is even quicker 
(~1 kyr) we do not expect many of the very massive stars to 
show up at a large radius (as red supergiants). Moreover, rapid 
rotation, which we do not take into account, may additionally 
limit the expansion of some of the most massive stars. 


3.2. Number of ULXs 


By the number of ULXs (#ULX), we understand the 
predicted number of observed systems at the present time. We 
find that #ULX strongly depends on the SFR _ history 
(Figure 2). 

Constant SFR naturally leads to a constant growth in #ULX, 
which starts early after the start of SF (~5 Myr) and after a few 
gigayears reaches a constant value of ##ULX = 9.3, 100, and 
56 for Zs, Z5/10, and Z;/100, respectively (see Section 4.2). 
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BHULXs appear first in both constant and burst SF 
scenarios. In the former case, the number of BHULXs increases 
to the age of a few hundred megayears, at which time it reaches 
a constant value of 3.8, 92, and 51 for Zo, Z5/10, and 
Ze /100, respectively. In the case of a burst SF, the formation 
time of BHULXs is limited nearly exactly to the duration of the 
burst (S200 Myr). 

NSULXs resulting from different evolutionary routes 
described in Section 3.1 appear sequentially in their chron- 
ological order and form the characteristic steep sections 
(constant SF) or humps (burst SF) at ~100 Myr (Ryns.ms), 
~200 Myr (Rwns,HG), and ~1 Gyr (Rns,rc). 

ULX populations younger than ~ 100 Myr are dominated by 
the BHULX. For constant SF, the NSULXs become more 
abundant than the BHULXs (58% of the ULX accretors) only if 
Z = Zs. For lower metallicities, the ##BHULX grows rapidly 
during the first few hundred megayears, whereas the rate of the 
NSULXs formation is much lower, and they are not able to 
match the abundance of the BHULXs during the next 10 Gyr. 
As a result, the NSULXs account for not more than ~12% and 
~8% of all ULXs for Z5/10 and Z,/100, respectively. 
However, for burst SF and all investigated metallicities, the 
NSULXs exceed the BHULXs in numbers within <1 Gyr, and 
after a few gigayears they become the only surviving ULXs. 


3.3. ULXs with Luminosities Exceeding 1 of erg gt 


Several surveys found a break in the X-ray luminosity 
function (XLF) at <2 x 10! erg s~!, which may point toward 
different populations below and above this luminosity (e.g., 
Grimm et al. 2003; Swartz et al. 2011; Mineo et al. 2012). 
Thus, it is interesting to point out that hyper-luminous X-ray 
sources (HLX) defined as ULXs with Lx = 104! erg s~!, thus 
clearly above the 2 x 10!°ergs~! threshold, have different 
donors than the ULXs with Lx < 10*! erg s~!. We find that 
52%-84% of the ULXs (described in Section 3) have 
luminosities in the range of (1-3) x 10°° ergs”. 

Approximately 90% of HLXs are BHULXs with HG donors. 
NS accretors in HLXs may obtain nearly as high luminosities 
as BH accretors and are typically accompanied by evolved 
helium stars (evHeS). It is noteworthy that the evolutionary 
HLX routes match those reported for HLXs_ with 
Lx > 10* erg s~! in Wiktorowicz et al. (2015). The current 
results, however, are based on a much wider range of models 
and a more physical treatment of NS/BH accretion, and we 
include here the HLX evolutionary paths for completeness 
(Table 2). Notice that the PULX NGC5907 ULX1 has a 
luminosity >104! erg s~' (Fiirst et al. 2016). 

15-60 Myr; Route Rgyuc: Typical companion HG—A 
typical system begins its evolution as a 25-50 M, primary with 
a 5-10 Ms companion. After ~5 Myr and the loss of 1-10 Me 
in stellar wind the primary fills the RL as a CHeB star and 
commences the CE phase. It loses a large fraction of its mass 
(~50%), but the separation shrinks to ~15 Ro. After additional 
~0.5 Myr a ~7-18 M. BH forms in a direct collapse. The 
secondary becomes an HG star ~10-55 Myr later and it 
expands rapidly. Shortly after, it fills its RL and starts the MT. 
For a short time (<0.1 Myr), the MT is high enough to power a 
ULX with Lx > 10+! erg s~! (see Figure 3, top). 

17-75 Myr; Route Rys,evues: Typical companion evHeS— 
The stars on ZAMS are less massive and have smaller mass 
ratios than in the case of Rgy yo. When the primary fills the RL 
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Table 4 
#ULX for Burst and Constant SFR 


Time Since the Start of Star Formation 


Metallicity 100 Myr 1 Gyr 
#ULX 4.0 x 10° 1.8 x 10! 
Fis #BHULX 3.7 x 107 2.8 x 1077 
#NSULX 3.5 x 10! 1.8 x 10! 
#ULX 73 x 10° 1.1 x 10° 
Zo/10 #BHULX 7.2 x 10° 4.0 x 10° 
#NSULX 8.1 x 10! 1.1 x 10° 
#ULX 5.0 x 10° 2.9 x 10! 
Z/100 ##BHULX 4.9 x 10° 
#NSULX 1.6 x 10! 2.9 x 10! 


Burst SFR" Constant SFR 

5 Gyr 10 Gyr 10 Gyr 
9.5 x 107! 1.1 x 107! 9.3 x 10° 
aoe any 3.8 x 10° 

9.5 x 107! 1.1 x 107! 5.5 x 10° 
7.0 x 107! 7.5 x 10-7 1.0 x 10° 
ae et 9.2 x 10! 

7.0 x 107! 7.5 x 10-7 1.3 x 10! 
1.0 x 107! 5.0 x 10-7 5.6 x 10! 
oe = 5.1 x 10! 

1.0 x 107! 5.0 x 107-7 5.0 x 10° 


Note. Number of ULXs (##ULX) in the reference model with a division on BHULXs (#4BHULX) and NSULXs (#NSULX) for different population ages. BHULX 
are also present late after the burst (~1 Gyr), but ##BHULX is negligible in that population’s age in comparison to ##4NSULX or the ##BHULX during the burst. 


* Burst duration is 100 Myr. 


Table 5 
#fULX(Lx min) 100 Myr after SF Start 
Lx min 
Metallicity 10° erg s~! 3 x 10° erg s—! 10“ erg s~! 104! erg s“! 
#ULX 4.0 x 10° 6.5 x 10! 1.4 x 10! 5.0 x 10°! 
Zs #BHULX 3.7 x 10° 4.0 x 10! 4.7 x 10° 3.8 x 10°! 
#NSULX 3.5 x 10! 2.5 x 10! 9.4 x 10° 13 x 107! 
#ULX 7.3 x 10° 1.9 x 103 2.0 x 10? 1.2 x 10! 
Zo [10 #BHULX 7.2 x 10° 1.9 x 10° 1.8 x 10° 1.1 x 10! 
#NSULX 8.1 x 10! 4.0 x 10! 1.2 x 10! 43 x 107! 
#ULX 5.0 x 10° 2.4 x 10° 3.9 x 10° 1.5 x 10! 
Z5/100 #BHULX 4.9 x 10° 2.4 x 10° 3.8 x 10° 1.5 x 10! 
#NSULX 1.6 x 10! 1.2 x 10! 9.0 x 10° 7.5 x 10° 


Note. Number of ULXs depending on the minimal luminosity (Lx min) for different metallicities with a division on BHULXs and NSULXs for ADI BK model. 


as an RG, the MT is stable and proceeds for a few tens of kyr. 
The donor losses most of its hydrogen envelope and after a few 
Myr at the age of ~25 Myr it forms an NS through an SN 
explosion. The companion, being an RG, expands and fills its 
RL. This time the MT is unstable due to a high-mass ratio and 
the CE occurs. Afterward, the secondary becomes a low-mass 
(a few Solar masses) helium star (HeS). The secondary 
continues its evolution for the next few Myr s, then expands 
and fills its RL again. This time, although being very strong, 
the MT is stable and for ~| kyr it powers an HLX (see 
Figure 3, bottom). 


4. Discussion 


ULXs are observed in all kinds of galaxies (e.g., Mush- 
otzky 2006). However, ULXs appear to be more abundant in 
the star-forming galaxies than in the elliptical galaxies (e.g., 
Irwin et al. 2004; Swartz et al. 2004; Feng & Soria 2011; Wang 
et al. 2016; however, see Soria 2007). In addition, ULXs are 
commonly detected in star-burst regions (e.g., Fabbiano et al. 
2001; Gao et al. 2003). Feng & Soria (2011) pointed out that 
this connection is particularly valid for the most luminous 


ULXs. On the contrary, Mushotzky (2006) performed a direct 
examination of the ULXs’ positions and found that approxi- 
mately one-third of them had to form outside of the star- 
forming regions. 

These observations may be understood on the ground of our 
results. On the one hand, the spatial correlation of ULXs and 
SF regions is expected because during the SF episodes ULXs 
form mainly through routes Rgyms and 7ns.ms, which 
produce numerous sources early after the start of SF. 
Additionally, the most luminous ULXs (HLXs) formed in 
our simulations evolve from heavier stars at an early age 
(tage S 200 Myr). On the other hand, the ULX phase in 
systems formed outside the SF regions occurs later on after 
ZAMS (Rns,yuc, 2ns.Rc)- However, we found that these late- 
time ULXs constitute only 1%-10% of all ULXs for 
constant SF. 


4.1. NS Accretors in ULXs 


Our results predict a large fraction of NS accretors among the 
ULXs. Current observational status is consistent with our results. 
The discovery of the first NSULX (Bachetti et al. 2014) proved 
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Table 6 
Number of ULXs and Maximal Luminosity (Constant SFR*) 


Model” Number per MWEG‘ Lx max (erg s~h4 
#ULX #BHULX #NSULX BHULX NSULX 
Zo 
ADO BN 3.1 x 10! 3.1 x 10° (10%) 2.8 x 10! (89%) 8.4 x 10 3.3 x 10% 
ADO BOl 1.8 x 10° 1.1 x 10° (0%) 1.8 x 107 (99%) 8.4 x 10° 3.3 x 10° 
ADO BK 2.2 x 10! 2.7 x 10° (12%) 1.9 x 10! (87%) 2.6 x 10% 1.0 x 10° 
ADO BS 2.0 x 10! 2.3 x 10° (11%) 1.8 x 10! (88%) 5.5 x 10% 2.2 x 10° 
AD1 BN 1.5 x 10! 4.4 x 10° (29%) 1.1 x 10' (70%) 3.0 x 10° 5.1 x 10° 
AD1 BO1 2.0 x 10? 1.5 x 10° (0%) 2.0 x 10° (99%) 3.0 x 104! 5.1 x 10% 
AD1 BK 9.3 x 10° 3.8 x 10° (41%) 5.5 x 10° (58%) 9.2 x 10” 1.6 x 10” 
AD1 BS 2.0 x 10! 2.9 x 10° (14%) 1.7 x 10! (85%) 2.0 x 10"! 3.4 x 10%° 
Ze/10 
ADO BN 1.6 x 10° 9.4 x 10! (58%) 6.7 x 10! (41%) 6.8 x 10% 2.3 x 10 
ADO BOl 3.7 x 10° 4.3 x 10' (11%) 3.3 x 10° (88%) 6.8 x 10% 2.3 x 10° 
ADO BK 1.3 x 10° 8.4 x 10! (63%) 5.0 x 10! (36%) 2.1 x 10% 7.1 x 10% 
ADO BS 1.1 x 10? 6.4 x 10! (60%) 4.1 x 10! (39%) 4.5 x 10% 1.5 x 10% 
AD1 BN 1.3 x 10° 1.1 x 107 (83%) 2.1 x 10! (16%) 8.6 x 10° 4.5 x 10° 
AD1 BOl 4.2 x 10° 5.0 x 10' (11%) 3.7 x 10° (88%) 8.6 x 104! 4.5 x 10% 
AD1 BK 1.0 x 10° 9.2 x 10! (87%) 1.3 x 10' (12%) 2.7 x 10% 1.4 x 10” 
AD1 BS 1.1 x 107 6.8 x 10' (61%) 4.2 x 10' (38%) 5.7 x 10" 3.0 x 10%° 
Ze /100 

ADO BN 8.5 x 10! 6.1 x 10! (72%) 2.3 x 10! (27%) 8.9 x 10% 2.5 x 10% 
ADO BOl 8.6 x 10! 1.4 x 10' (16%) 7.2 x 10! (83%) 8.9 x 10% 2.5 x 10° 
ADO BK 7.2 x 10! 5.6 x 10! (78%) 1.6 x 10! (21%) 2.7 x 10% 7.7 x 10% 
ADO BS 4.9 x 10! 3.8 x 10! (78%) 1.1 x 10! (21%) 5.9 x 10° 1.6 x 10° 
AD1 BN 6.6 x 10! 5.9 x 10! (89%) 7.2 x 10° (10%) 1.3 x 107 5.5 x 10°? 
AD1 BOl 8.0 x 10! 1.6 x 10' (19%) 6.4 x 10! (80%) 1.3 x 10” 5.5 x 10% 
AD1 BK 5.6 x 10! 5.1 x 10! (91%) 5.0 x 10° (8%) 4.0 x 10% 1.7 x 10” 
AD1 BS 5.0 x 10! 3.7 x 10! (73%) 1.3 x 10! (26%) 8.6 x 104! 3.6 x 10% 
Notes. 


* At the age of 10 Gyr. 


> ADO /1—“upper limit’ /logarythmic accretion model (Appendix A.1 and Section 2.1); BN/B0O1/BK/BS—beaming models (Appendix A.2 and Section 2.2). 


© Milky Way Equivalent Galaxy. 


4 The maximal obtained luminosity for the particular accretor. In all cases, the beaming that corresponds to these luminosities is saturated, so the beaming factor 
b =1,0.1, ~3.2 x 10-3, and ~0.15 for BN, B01, BK, and BS, respectively. Corresponding opening angles are 180°, ~52°, ~9°, and ~64°, respectively. 


that apparent luminosities orders of magnitude larger than the EL 
are possible in nature. Kluzniak & Lasota (2015) proposed that 
the M82 X-1 and the P13 in NGC 7793 may also be PULXs. 
The following discoveries (Fiirst et al. 2016; Israel et al. 2017a) 
showed that NSs may indeed be widespread among the ULXs. 
Finally, a known X-ray pulsar SMC X-3 was reported to reach 
an outburst luminosity of 2.5 x 10°° erg s~! (Tsygankov et al. 
2017). Pintore et al. (2017) argued that the energy spectra of 
several ULXs can be explained with a “pulsator-like” continuum 
model used to describe Galactic XRBs containing NSs (e.g., 
White et al. 1995; Coburn et al. 2001). 

Theoretical considerations further support our findings. King 
& Lasota (2016) showed that a ULX with a weakly magnetized 
NS and a ULX with a BH may be observationally hard to 
distinguish if the NS does not pulsate. King et al. (2017) argued 
that the conditions for pulsations in PULXs are relatively strict 
and a large population of NSULXs will mimic the BHULXs. 

The significant number of NSULX that follows from most of 
our simulated models may be understood on the statistical 
ground. While it is easier to obtain large MT rates in the 
BHULXs, the NSULXs progenitors are more abundant in the 
ZAMS distributions. In general, NSs form from stars with 


ZAMS masses ~8-22 Mz, whereas BHs form from more 
massive stars with ~22-150 M,. For aime = —2.7, it gives the 
ratio of NS to BH progenitors of ~4.8, and for ayyr = —2.3 
this ratio is 3.0. Moreover, in the NSULX route, which 
dominates in the old stellar populations (Rns.rG), the primary 
on ZAMS may have an even lower mass of ~6 Ms, which 
gives the ratio of NS to BH progenitors of 4.8 for 
Qimr = —2.3. It first becomes a WD before accreting mass 
and collapsing into an NS. Additionally, donors in NSULXs 
are less massive than in BHULXs and, therefore, also more 
abundant on ZAMS. We note that low-mass donors may occur 
also in BHULXs, but it is far harder (and therefore less 
probable) for them to survive the CE. Even if they succeeded, 
the post-CE separation is so small that the RLOF occurs earlier 
than what is typical for NSULXs. All of these results are for a 
significantly smaller number of BHULXs with low-mass 
(S2 Mz) donors than in NSULXs. 

During short SF episodes, the ULX population is dominated 
by sources with BH accretors. BHs may acquire stable MT 
from heavier companions than NSs. Heavier stars expand more 
and faster, therefore, BHULXs age during ULX phase is of the 
order of ~10 Myr, whereas for NSULX it is ~100 Myr (route 
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Table 7 
General Parameters of the Simulational Results (Burst SFR 100 Myr Ago) 


Number per MWEG* 


Lx max (erg Ss” ye 


Model* 

#ULX #BHULX #NSULX BHULX NSULX 

Zo 
ADO BN 3.4 x 10? 2.8 x 10° (83%) 5.8 x 10! (16%) 8.4 x 10“ 3.3 x 10% 
ADO BOl 4.7 x 10° 9.3 x 10' (19%) 3.8 x 10° (80%) 8.4 x 10° 3.3 x 10° 
ADO BK 2.8 x 10° 2.6 x 107 (90%) 2.8 x 10! (9%) 2.6 x 10" 1.0 x 10% 
ADO BS 2.8 x 10° 2.1 x 10° (75%) 6.8 x 10! (24%) 5.5 x 10% 2.2 x 10° 
AD1 BOl 5.7 x 10° 1.3 x 10% (22%) 4.4 x 10° (77%) 3.0 x 104! 5.1 x 10*° 
AD1 BK 4.0 x 10° 3.7 x 107 (91%) 3.5 x 10! (8%) 9.2 x 10” 1.6 x 10” 
AD1 BS 3.8 x 10° 2.8 x 10° (75%) 9.4 x 10! (24%) 2.0 x 1071 3.4 x 10% 
Zo/10 
ADO BN 7.5 x 10° 7.2 x 10° (96%) 2.9 x 10? (3%) 6.8 x 10% 1.2 x 10% 
ADO BOl 3.0 x 10° 3.0 x 10° (98%) 4.9 x 10! (1%) 6.8 x 10% 1.2 x 10° 
ADO BK 6.6 x 10° 6.4 x 10° (98%) 1.2 x 10° (1%) 2.1 x 10% 3.7 x 10% 
ADO BS 5.2 x 10° 5.1 x 10° (98%) 6.7 x 10! (1%) 4.5 x 10% 7.9 x 10“ 
AD1 BN 8.6 x 10° 8.5 x 10° (98%) 1.0 x 107 (1%) 8.6 x 10°° 4.5 x 10°° 
AD1 BOl 3.9 x 10° 3.8 x 10° (96%) 1.4 x 107 (3%) 8.6 x 104! 4.5 x 10% 
AD1 BK 7.3 x 10° 7.2 x 10° (98%) 8.1 x 10' (1%) 2.7 x 10% 14 x 10” 
AD1 BS 5.7 x 10° 5.6 x 10° (97%) 1.5 x 10° (2%) 5.7 x 10"! 3.0 x 10%° 
Zs /100 

ADO BN 5.9 x 10° 5.9 x 10° (99%) 2.8 x 10! (0%) 8.9 x 10° 2.1 x 10% 
ADO BOl 13 x 10° 1.3 x 10° (99%) 4.4 x 10° (0%) 8.9 x 10*° 2.1 x 104% 
ADO BK 5.4 x 10° 5.4 x 10° (99%) 1.1 x 10! (0%) 2.7 x 10% 6.5 x 10° 
ADO BS 3.7 x 10° 3.7 x 10° (99%) 6.9 x 10° (0%) 5.9 x 10° 14 x 10“ 
AD1 BN 5.7 x 10° 5.7 x 10° (99%) 1.4 x 10! (0%) 1.3 x 10"! 5.5 x 10° 
AD1 BOl 1.5 x 10° 1.5 x 10° (96%) 4.9 x 10! (3%) 1.3 x 10% 5.5 x 10% 
ADI BK 5.0 x 10° 4.9 x 10° (99%) 1.6 x 10! (0%) 4.0 x 10% 1.7 x 10” 
AD1 BS 3.7 x 10° 3.6 x 10° (98%) 5.8 x 10! (1%) 8.6 x 104! 3.6 x 10%° 


Note. The same as in Table 6, but for burst SF, which started 100 Myr ago and lasted for 100 Myr. 


Rwns.ms)- AS a result BHULXs appear in larger numbers in the 
first ~ 100 Myr. The ULX phase may appear in NSULX as 
early as ~10 Myr after ZAMS. However, it demands a very 
specific pre-CE conditions, which occur in a relatively narrow 
ZAMS parameter space. For an average post-CE separation in 
early NSULXs (Rys_ms) the star needs ~ 100 Myr to fill its RL. 

The picture is different for the constant SF or long bursts. As 
the SF continues, the ##BHULX saturates for ~200 Myr after 
the SF start, because at this time a comparable number of 
sources begin and end their ULX phase. However, a significant 
number of NSULXs have their ULX phases during the post- 
MS expansion of their light donors. Therefore, the ##+NSULX 
grows continuously for several gigayears and new NSULX 
formation routes appear successively (see Section 3.1 for the 
sequence). For Z = Z,, the ##NSULX becomes higher than 
#tBHULX after ~1 Gyr. However, for lower metallicities, 
##NSULX does not overcome ##BHULX within 10 Gyr. 

When the SF is extinguished, gyms quickly disappear as 
their massive donors quickly end their lives (C500 Myr). On 
the other hand, NSULXs with low-mass donors may 
commence the ULX phase after several gigayears of evolution 
when the secondary fills the RL during the post-MS expansion 
(Rns.uG and Rys.rG). Therefore, after the end of the SF, 
NSULXs quickly (in ~100 Myr) become the dominant ULXs 
and after ~1 Gyr BHULXs disappear completely. 

The population synthesis of NSULXs was recently per- 
formed with the use of the BSE code by Fragos et al. (2015). 
However, they omitted the pre-SN evolution and explored a far 
smaller parameter space than that presented in this work. Even 


though they utilized the MESA code to calculate the MT rates, 
their maximum values are similar to ours (M < 10-7 Myr). 
We obtained less massive companions in NSULXs, but (Fragos 
et al. 2015) focused on M82 X-2 source with companion mass 
25.2 Ms (Bachetti et al. 2014) and assumed that heavier stars 
are necessary to survive the CE phase. 

Shao & Li (2015) also took the population synthesis 
approach to study NSULXs and found a significant population 
of these sources. They adopted a realistic envelope’s binding 
energy from Xu & Li (2010), but considered only a solar 
metallicity and constant beaming of b = 0.1. 


4.2. Dependence of #ULX on Metallicity 


Observations show that ULXs are often associated with low 
metallicity environments (e.g., Pakull & Mirioni 2002; Soria 
et al. 2005; Luangtip et al. 2015; Mapelli et al. 2010). Many 
authors have demonstrated that in a low metallicity environ- 
ment it is easier to form a massive BH (e.g., Mapelli et al. 
2009; Zampieri & Roberts 2009; Belczynski et al. 2010), and it 
is more common to obtain an RLOF onto a compact object 
(Linden et al. 2010). 

In this study, we find that the metallicity has a strong impact 
on #ULX, but only in SF regions (See Table 4). We show that 
the number of ##BHULX increases significantly for Z,/10 in 
comparison to Zs. However, the #NSULX is virtually 
unaffected by metallicity. Interestingly, for models with the 
lowest investigated metallicity (Z,/100), we found fewer 
ULXs than for Z = Z,/10. Therefore, our results suggest that 
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Table 8 
General Parameters of the Simulational Results (Burst SFR 1 Gyr Ago) 


Number per MWEG* 


Lx max (erg she 


Model* 

#ULX #BHULX #NSULX BHULX NSULX 

Zo 
ADO BN 8.3 x 10! 8.5 x 1077 (0%) 8.3 x 10' (99%) 2.0 x 10 2.2 x 10% 
ADO BOl 3.9 x 107 3.6 x 10° | (0%) 3.9 x 10? (99%) 2.0 x 10% 2.2 x 10% 
ADO BK 6.0 x 10! 3.5 x 10 + (0%) 6.0 x 10! (99%) 6.2 x 10“ 6.8 x 10 
ADO BS 5.5 x 10! 1.5 x 10-7 (0%) 5.5 x 10! (99%) 1.3 x 10% 14 x 10% 
ADI BN 3.8 x 10! 8.7 x 10° (0%) 3.8 x 10! (99%) 3.7 x 10°° 3.9 x 10° 
ADI BOl 5.2 x 10° 1.2 x 10°' (0%) 5.2 x 10° (99%) 3.7 x 10° 3.9 x 10% 
ADI BK 1.8 x 10! 2.8 x 10-7 (0%) 1.8 x 10' (99%) 1.1 x 107 1.2 x 10” 
ADI BS 8.0 x 10! 1.3 x 10° (0%) 8.0 x 10! (99%) 2.4 x 10% 2.6 x 10° 
Zo /10 
ADO BN 6.8 x 10° 4.1 x 10° (0%) 6.8 x 10° (99%) 4.8 x 10” 1.4 x 10% 
ADO BOl 8.5 x 107 1.2 x 10' (1%) 8.4 x 10 (98%) 48 x 10% 1.4 x 10“ 
ADO BK 5.7 x 10° 3.9 x 10° (0%) 5.7 x 10° (99%) 1.5 x 10° 43 x 10° 
ADO BS 2.8 x 10° 4.2 x 10° (1%) 2.8 x 10? (98%) 3.2 x 10% 9.2 x 10% 
ADI BN 1.2 x 10° 4.0 x 10° (3%) 1.2 x 10° (96%) 9.0 x 10° 3.5 x 10° 
ADI BOl 9.9 x 10° 5.0 x 10° (0%) 9.8 x 10° (99%) 9.0 x 10° 3.5 x 10% 
ADI BK 1.1 x 10° 4.0 x 10° 3%) 1.1 x 10° (96%) 2.8 x 10” 1.1 x 10” 
ADI BS 2.9 x 10° 3.8 x 10° (1%) 2.9 x 10° (98%) 5.9 x 10% 2.3 x 10% 
Zo /100 

ADO BN 2.2 x 10? 5.6 x 10-7 (0%) 2.2 x 10? (99%) 3.4 x 10 7.3 x 10% 
ADO BOl 1.7 x 10° 1.6 x 10° (0%) 1.7 x 10° (99%) 3.4 x 10% 73 x 10% 
ADO BK 1.6 x 10° 1.3 x 10-7 (0%) 1.6 x 10° (99%) 1.0 x 10% 2.2 x 10° 
ADO BS 7.2 x 10! 1.1 x 10-7 (0%) 7.2 x 10! (99%) 2.2 x 108 4.8 x 10% 
AD1 BN 3.1 x 10! 0.0 x 10° (0%) 3.1 x 10! (100%) 0.0 x 10° 2.8 x 10°? 
ADI BOl 2.2 x 10° 0.0 x 10° (0%) 2.2 x 10° (100%) 0.0 x 10° 2.8 x 10% 
AD1 BK 2.9 x 10! 0.0 x 10° (0%) 2.9 x 10! (100%) 0.0 x 10° 8.6 x 104 
ADI BS 6.3 x 10! 0.0 x 10° (0%) 6.3 x 10! (100%) 0.0 x 10° 1.8 x 10% 


Note. The same as in Table 6, but for burst SF, which started 1 Gyr ago and lasted for 100 Myr. 


the relation ##ULX(Z) is not monotonic (see Prestwich et al. 
2013). We plan to study this outcome in more detail in a 
forthcoming publication (G. Wiktorowicz et al. 2017, in 
preparation). 


4.3. Counterparts 


Several searches were performed to directly observe the 
ULX companions. However, only a handful have been detected 
among hundreds of known ULXs. Additionally, ULXs, being 
extra-galactic objects, are observationally biased toward more 
luminous and therefore more massive companions. In our 
results, NS counterparts in ULX are predominantly low-mass 
(S1.5 Mz) stars. If located in other galaxies, such objects 
would be difficult to detect. This may explain why it is hard to 
find counterparts to the majority of the ULXs. 

Some ULXs were observed near OB stars (e.g., Soria et al. 
2005). Such stars are too massive to provide a stable MT to an 
NS accretor, so they may only be accompanied by BH 
accretors, unless the MT occurs through a stellar wind (not 
considered in our study). OB stars evolve quickly (< 100 Myr), 
which suggests that they should be spatially associated with the 
SF regions, where the BHULXs are more numerous than the 
NSULXs according to our results. Typical BHULX compa- 
nions in our simulations are MS stars of 5.6-11 Mz, which 
matches the B spectral type. 

Roberts et al. (2008) found five potential optical ULX 
companions with the use of the Hubble Space Telescope. The 


optical search for companions in ULXs was performed also by 
Gladstone et al. (2013) who used the same instrument and 
found 13 + 5 counterparts to ULXs located within a distance 
of D < 5 Mpc. Most of them were MS stars or RGs. A near- 
infrared search was performed by Heida et al. (2014, 2016) 
who observed 62 ULXs in 37 galaxies located within 10 Mpc. 
They found 17 candidate counterparts, 13 of which are red 
supergiants. We note that red supergiants are present in our 
results in small numbers (< 1%). The relatively large number of 
these stars in observations may be an observational bias, as 
they are significantly brighter than the most typical companions 
predicted by our study (<2 M,). Nevertheless, red supergiants 
may be an interesting topic for a separate study. 

The companion to NSULX P13 in NGC 7793 was detected 
and estimated to be a BI9a star with mass 18-23 M, (Motch 
et al. 2014). The radius of such a star is 4-700 Ro, depending 
on its evolutionary stage. An orbital period estimate 
(Pxb = 64 days) puts an RL radius of a donor at ~120 Ro 
(under the assumption of circular orbit and Mys = 1.4 Mz). 
The BI9a star may obtain such a radius during the CHeB phase. 
According to a standard prescription, in a high mass-ratio 
system, the RLOF will be unstable and the CE will occur. 
However, Pavlovskii et al. (2017) found that the range of stable 
mass ratios may be wider than previously thought. We will 
analyze this possibility in a separate study (Wiktorowicz et al. 
2017, in preparation). BHULXs with high-mass donors may be 
predecessors of double BH mergers (Belczynski et al. 2016). 
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Table 9 
General Parameters of the Simulational Results (Burst SFR 5 Gyr Ago) 


Number per MWEG* 


Lx, max (erg Ss” ae 


Model* 

#ULX #BHULX #NSULX BHULX NSULX 

Zo 
ADO BN 7.3 x 10° 0.0 x 10° (0%) 7.3 x 10° (100%) 0.0 x 10° 7.2 x 10% 
ADO BOl 6.6 x 10! 0.0 x 10° (0%) 6.6 x 10! (100%) 0.0 x 10° 7.2 x 10% 
ADO BK 4.9 x 10° 0.0 x 10° (0%) 4.9 x 10° (100%) 0.0 x 10° 2.2 x 10% 
ADO BS 3.5 x 10° 0.0 x 10° (0%) 3.5 x 10° (100%) 0.0 x 10° 4.7 x 10“ 
ADI BN 72x 107! 0.0 x 10° (0%) 7.2 x 10-' (100%) 0.0 x 10° 2.9 x 10° 
AD1 BO1 5.0 x 10! 0.0 x 10° (0%) 5.0 x 10! (100%) 0.0 x 10° 2.9 x 10% 
ADI BK 9.5 x 107! 0.0 x 10° (0%) 9.5 x 10°! (100%) 0.0 x 10° 8.9 x 104! 
AD1 BS 2.0 x 10° 0.0 x 10° (0%) 2.0 x 10° (100%) 0.0 x 10° 1.9 x 10” 
Zo/10 
ADO BN 6.7 x 10° 0.0 x 10° (0%) 6.7 x 10° (100%) 0.0 x 10° 9.4 x 10” 
ADO BOl 6.5 x 10! 7.2 x 10°' (1%) 6.4 x 10! (98%) 6.1 x 10°° 9.4 x 10% 
ADO BK 4.9 x 10° 0.0 x 10° (0%) 4.9 x 10° (100%) 0.0 x 10° 2.9 x 10% 
ADO BS 44 x 10° 0.0 x 10° (0%) 4.4 x 10° (100%) 0.0 x 10° 6.2 x 10% 
ADI BN 1.4 x 10° 0.0 x 10° (0%) 1.4 x 10° (100%) 0.0 x 10° 3.0 x 10°? 
ADI BOl 5.1 x 10! 0.0 x 10° (0%) 5.1 x 10! (100%) 0.0 x 10° 3.0 x 10% 
AD1 BK 7.0 x 107! 0.0 x 10° (0%) 7.0 x 107! (99%) 0.0 x 10° 9.2 x 104! 
ADI BS 3.2 x 10° 0.0 x 10° (0%) 3.2 x 10° (100%) 0.0 x 10° 2.0 x 10*° 
Ze /100 

ADO BN 1.9 x 10° 0.0 x 10° (0%) 1.9 x 10° (100%) 0.0 x 10° 1.2 x 10% 
ADO BOl 1.9 x 10! 2.0 x 10°-* (0%) 1.9 x 10! (99%) 6.8 x 10°? 1.2 x 10“ 
ADO BK 1.4 x 10° 0.0 x 10° (0%) 1.4 x 10° (100%) 0.0 x 10° 3.7 x 10° 
ADO BS 2.0 x 10° 0.0 x 10° (0%) 2.0 x 10° (100%) 0.0 x 10° 7.9 x 10% 
ADI BN 3.0 x 10-4 0.0 x 10° (0%) 3.0 x 10-4 (100%) 0.0 x 10° 3.0 x 10°? 
ADI BOl 4.9 x 10° 1.2 x 10°' (2%) 4.8 x 10° (97%) 6.5 x 10°? 3.0 x 10% 
ADI BK 1.0 x 107! 0.0 x 10° (0%) 1.0 x 10°! (100%) 0.0 x 10° 9.2 x 1071 
ADI BS 7.9 x 107! 0.0 x 10° (0%) 7.9 x 10~' (100%) 0.0 x 10° 2.0 x 10% 


Note. The same as in Table 6, but for burst SF which started 5 Gyr ago and lasted for 100 Myr. 


5. Summary and Conclusions 


We analyzed accretion onto compact stars in order to investigate 
the parameters and characteristics of the ULX population. We 
demonstrate that on the grounds of the current understanding of the 
accretion physics and stellar evolution we are able to reproduce the 
vast majority of the observed population of ULXs without the 
need for intermediate-mass black hole accretors. 

We performed a population synthesis study with the use of 
the StarTrack population synthesis code to simulate 
2 x 10’ isolated binaries with initial parameters leading to 
the formation of XRBs (both LMXB and HMXB). Our 
simulation grid involves various accretion models, beaming 
prescriptions, and metallicities. 

Our main conclusions are as follows. 


(1) ULX with NS accretors dominate the post-burst ULX 
populations, and constant SF (duration > 1 Gyr) in high-Z 
environments, which is a natural consequence of the 
current understanding of the binary evolution. NSULXs 
are present in significant numbers (= 10%) also during the 
SF bursts and in lower-Z ULX populations. 

(2) ULXs appear in a very specific sequence after the start of 
the SF (t depicts here the time since the beginning of 
the SF): 

1. t= 4-40 Myr BH-MS (5.6-11 M3), 

2. t © 6-800 Myr NS—MS (0.9-1.5 Mo), 

3. t & 430-1100 Myr NS-HG (0.6-1.0 M,), 
4. t = 540-4400 Myr NS-RG (~1.0 M,). 
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(3) We found that NSULXs may reach luminosities as high as 
those of BHULXs (Lx.max > 104! erg s~!). 

(4) The most luminous ULXs (Lx = 10*! erg s~!) contain HG 
donors (BHULXs; My, & 1.2-3.7 Mas) or evolved helium 
stars (NSULXs; M, = 1.7-2.6 Ms), which overfill their 
RL and transfer mass on thermal-timescale. They form 
typically within 15—75 Myr since ZAMS. 

(5) The number of ULXs is anti-correlated with metallicity. 
However, for very low-Z this relation changes sign. 
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ported in part by a grant from the French National Space Center 
CNES. K.B. acknowledges support from the NCN grants: 
OPUS 2015/19/B/ST9/03188 and MAESTRO 2015/18/A/ 
ST9/00746. G.W. was supported by the NCN grant OPUS 
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Table 10 
General Parameters of the Simulational Results (Burst SFR 10 Gyr Ago) 
Model" Number per MWEG* Lx. max (erg s~)? 
#ULX #BHULX #NSULX BHULX NSULX 
Zo 
ADO BN 1.2 x 10° 0.0 x 10° (0%) 1.2 x 10° (100%) 0.0 x 10° 8.7 x 107 
ADO BO1 8.1 x 10° 0.0 x 10° (0%) 8.1 x 10° (100%) 0.0 x 10° 8.7 x 10% 
ADO BK 8.3 x 107! 0.0 x 10° (0%) 8.3 x 107! (100%) 0.0 x 10° 2.7 x 10% 
ADO BS 74 x 107! 0.0 x 10° (0%) 74 x 10~' (100%) 0.0 x 10° 5.7 x 10% 
AD1 BN 2.9 x 10°* 0.0 x 10° (0%) 2.9 x 10°* (100%) 0.0 x 10° 3.0 x 10° 
AD1 BO1 2.6 x 10° 0.0 x 10° (0%) 2.6 x 10° (100%) 0.0 x 10° 3.0 x 10%° 
AD1 BK 1.1 x 107! 0.0 x 10° (0%) 1.1 x 107! (100%) 0.0 x 10° 9.2 x 104! 
AD1 BS 3.6 x 10°! 0.0 x 10° (0%) 3.6 x 10°' (100%) 0.0 x 10° 2.0 x 10% 
Z5/10 
ADO BN 1.5 x 10° 0.0 x 10° (0%) 1.5 x 10° (100%) 0.0 x 10° 5.4 x 10” 
ADO BO1 1.2 x 10! 0.0 x 10° (0%) 1.2 x 10! (100%) 0.0 x 10° 54x 10% 
ADO BK 1.1 x 10° 0.0 x 10° (0%) 1.1 x 10° (100%) 0.0 x 10° 1.7 x 10° 
ADO BS 9.4 x 107! 0.0 x 10° (0%) 9.4 x 107! (100%) 0.0 x 10° 3.6 x 108 
AD1 BN 2.1 x 107! 0.0 x 10° (0%) 2.1 x 107! (100%) 0.0 x 10° 2.8 x 10°? 
AD1 BO1 7.3 x 10° 0.0 x 10° (0%) 7.3 x 10° (100%) 0.0 x 10° 2.8 x 10° 
AD1 BK 7.5 x 10° 0.0 x 10° (0%) 7.5 x 10°* (100%) 0.0 x 10° 8.6 x 107 
ADI BS 2.9 x 107! 0.0 x 10° (0%) 2.9 x 107! (100%) 0.0 x 10° 1.8 x 10” 
Zs /100 

ADO BN 6.2 x 107! 0.0 x 10° (0%) 6.2 x 107! (100%) 0.0 x 10° 5.1 x 10” 
ADO BOl 5.9 x 10° 0.0 x 10° (0%) 5.9 x 10° (100%) 0.0 x 10° 5.1 x 10% 
ADO BK 4.4 x 107! 0.0 x 10° (0%) 4.4 x 107! (100%) 0.0 x 10° 1.6 x 10° 
ADO BS 3.8 x 107! 0.0 x 10° (0%) 3.8 x 10-!| (100%) 0.0 x 10° 3.4 x 10% 
AD1 BN 1.5 x 107+ 0.0 x 10° (0%) 1.5 x 10-4 (100%) 0.0 x 10° 2.8 x 10°? 
AD1 BO1 2.9 x 10° 0.0 x 10° (0%) 2.9 x 10° (100%) 0.0 x 10° 2.8 x 10%° 
ADI BK 5.0 x 1077 0.0 x 10° (0%) 5.0 x 1077 (100%) 0.0 x 10° 8.6 x 107! 
AD1 BS 2.3 x 10> 0.0 x 10° (0%) 2.3 x 10> (100%) 0.0 x 10° 1.8 x 10% 


Note. The same as in Table 6, but for burst SF, which started 10 Gyr ago and lasted for 100 Myr. A total absence of BHULXs may be observed. 


0707. This research was supported in part by the National 
Science Foundation under grant No. NSF PHY-1125915. 


Appendix 
Details of the Simulations 


H/R 


If not stated differently, throughout the article a “luminosity” 
means an X-ray luminosity in the 0.2-10 keV band. We do not 
take into account NS magnetic fields and include only non- 
rotating BHs. All of the frequently used abbreviations are 
summarized in Table 3. 

We utilized the EL for luminosity (Zgaa) and mass accretion 
( Mga) as or 1 10 100 “4000 

. MI Meda 
(4s), (4) Figure 4. Ratio of photosphere height to disk radius (H/R) in relation to the 

Mos mass transfer rate (in Eddington units). Orange squares represent the results of 


GRRMHD simulations with the use of KORAL code (A. Sadowski et al. 2016, 
private communication). The blue line shows the fit to the data (Equation (8)). 


Lead = 9 Mpaac? = 2.51 x 1078 erg s~! 
Edd = 7) Medd g 14x 


where 7 = 1/12 for BH accretors and ~0.2 for NSs, X is a 


fraction of hydrogen in the envelope of the donor. X = 0.7 for A.1. The “Upper Limit” Accretion Model (ADO) 
hydrogen reach donors and X = 0 for HeS and WD. . : : ; : 

Most of the results are provided as present time distributions. In addition to the accretion model described in Section 2.1, 
This means the distributions as they will be visible by an to which we will refer to as AD1, we investigated a possibility 
observer in present time (10 Gyr after the SF beginning for that all mass may be transferred from the donor and efficiently 
constant SF, or (100 Myr, | Gyr, 5 Gyr, and 10 Gyr after the SF accreated by the accretor, 1.e., 
beginning for burst SF). We omit issues related to obscuration Mess = Mabe. (5) 
or light propagation time. 

Below, we list the additional models tested in our grid of This corresponds to the highest accretion rate that could 
simulations. The standard model will be referred to as AD1 BK. potentially occur in a given system, but is most likely not 
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Figure 5. Evolution of the relative abundance of young (fage < 50 Myr; mainly routes Rgy.ms and 7ns.ms) and old (fage > 500 Myr; mainly routes Rys,Ho and 
Rns,rc) ULXs through time for constant star-formation rate and two metallicities: Z = Ze (left) and Z = Z,/10 (right). 


realistic for high Mryor. Therefore, this model should be 
considered as a rough upper limit for accretion. 
The X-ray luminosity is calculated as 


€G Mpuyns Mrior 

Race 
where € is a conversion efficiency of gravitational energy into 
radiation equal 1.0 for an NS (surface accretion) and 0.5 for a 
BH (disk accretion), Mgyyns is the mass of the accretor, Race 
is the radius of an NS (assumed to be 10 km) or a BH (3 
Schwarzschild radii, non-rotating BH), and 7 is the radiative 
efficiency of a standard thin disk equal ~0.2 for an NS and 1/ 
12 for a BH (e.g., Shakura & Sunyaev 1973). 


ig= (6) 


= y 2 
7] Mriorc’, 


A.2. Models of the Beaming 


In addition to the beaming model described in Section 2.2 to 
which we will refer to as BK, we investigated also the 
following models: 


A.2.1. No Beaming (BN) 
For this model, we assumed isotropic emission, therefore, 
b=1 (7) 


for all systems and all MT rates. 


A.2.2. Constant Beaming, no Saturation (BO/) 


We start with the simplest model of constant beaming, and 
we apply it to all sources. We consider only one case of 
b = 0.1 (6 & 52°) to compare it with more realistic prescrip- 
tions described below. Such a constant beaming will lower the 
total luminosity, which is required for a source to be observed 
as a ULX to b x 10°’ erg s~!. This will generally increase the 
number of predicted ULX sources. On the other hand, Rps,b 
will be lower for all systems, which will decrease this number. 
These two processes may make the predicted number of ULXs 
lower or higher in comparison to models without beaming. 


A.2.3. Sqdowski’s Model with Saturation (BS) 


This beaming prescription is based mainly on the results of 
theoretical and numerical analyses presented in Lasota et al. 
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(2016) and based on the results of KORAL (Sadowski et al. 
2014) GRRMHD simulations. They showed that super- 
Eddington disks never become geometrically thick, because 
the thickness of the slim disk does not depend on the MT rate. 
Even for very high si the ratio of photosphere height (H) to 
radius (R) is H/R < 1.6. We fitted a phenomenological model 
to approximate the relation between H/R and m found by 
Lasota et al. (2016; Figure 4), and we obtained 


1.6 (8) 


for non-rotating BH accretion with non-saturated magnetic 
field. This relation corresponds to R = 30 Rg, but H/R should 
not be significantly larger for other radii. The equation shows a 
moderate photosphere height also for sub-Eddington MT rates. 
We used the same prescription for NS accretors. 

The H/R is related to the opening angle 0 as 


Hy 0 
= = tan—, 9 
(E) = 95 ° 
which allows us to derive the beaming factor as 
0 
b=1-—cos-, 10 
5 (10) 


and the apparent isotropic luminosity is calculated as in 
Equation (2). 

The beaming becomes nearly constant for m = 100, which 
gives the minimum opening angle Onin © 64° (b © 0.15). The 
prescription works for the entire range of 7 and shows a small 
collimation also for nearly Eddington MT (m < 1). However, 
below m#0.12 the beaming factor is always greater 
than b > 0.95(0 > 170°). 


A.2.4, Additional Characteristics of the Simulated Populations 
of ULXs 


Table 5 contains the dependence of #ULX on Lyx min. The 
distributions of the general population’s parameters are in 
Tables 6-10. The evolution of the ##ULX through population 
age are in Figures 6 and 7. Tables 11-15 contain the 
distributions of companions. 
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Z,,/100 AD1 BK 


2 3 a al 
0 10 10° 10 
t [Myr] 
Figure 6. Relation between the number of ULXs and the age of the population for constant SFR of 6.0 Ms yr! (total stellar mass equal Milky Way Equivalent 


Galaxy). The thick line corresponds to BHULXs, dashed to NSULXs, and solid to all ULXs. Ze depicts the solar metallicity, ADO/1 are the accretion models 
(Appendix A.1 and Section 2.1) and BN, BO1, BK, and BS stand for beaming models (Appendix A.2 and Section 2.2). 
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Figure 7. Same as in Figure 6, but for burst star formation with a duration of 100 Myr and SFR of 600 Ms yr~’. 
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Z,/100 ADO B1 


Z,/100 ADO BO1 


Z,,/100 ADO BK 


Z,,/100 ADO BS 


Z,,/100 AD1 B1 


Z,,/100 AD1 BOL 


Z,,/100 AD1 BK 


Z,,/100 AD1 BS 


1 


ST 


Model 


ADO BN 
ADO BO1 
ADO BK 
ADO BS 
AD1 BN 
AD1 BOL 
AD1 BK 
AD1 BS 


ADO BN 
ADO BOL 
ADO BK 
ADO BS 
AD1 BN 
AD1 BOL 
AD1 BK 
AD1 BS 


ADO BN 
ADO BO1 
ADO BK 
ADO BS 
AD1 BN 
AD1 BOL 
AD1 BK 
AD1 BS 


MS 


0.82 
0.91 
0.91 
0.92 
0.85 
0.94 
0.94 
0.94 


0.68 
0.87 
0.74 
0.79 
0.69 
0.91 
0.76 
0.81 


0.77 
0.89 
0.83 
0.87 
0.76 
0.91 
0.84 
0.89 


HG 


0.18 
0.07 
0.09 
0.08 
0.15 
0.05 
0.06 
0.06 


0.30 
0.10 
0.24 
0.20 
0.28 
0.08 
0.21 
0.17 


0.09 
0.04 
0.05 
0.05 
0.10 
0.04 
0.06 
0.05 


RG 


0.01 
0.02 


0.01 
0.01 


0.01 
0.01 


CHeB 


0.01 


0.01 
0.01 
0.02 


0.01 
0.01 


0.13 
0.06 
0.12 
0.08 
0.13 
0.05 
0.10 
0.06 


HeMS 


0.01 


0.01 


HewWD 


Table 11 


Companions (Constant SFR) 


HybWD 


MS 


0.03 
0.02 
0.01 
0.01 
0.04 
0.05 
0.13 
0.09 


0.06 
0.02 
0.03 
0.03 
0.10 
0.06 
0.11 
0.08 


0.02 
0.03 
0.01 
0.02 
0.04 
0.13 
0.12 
0.10 


HG 


Zo 
0.05 
0.02 
0.04 
0.06 
0.10 
0.04 
0.10 
0.16 

Zo/10 
0.47 
0.03 
0.51 
0.49 
0.39 
0.04 
0.54 
0.52 

Zo/100 
0.48 
0.03 
0.48 
0.29 
0.45 
0.05 
0.47 
0.36 


RG 


0.05 
0.01 
0.04 
0.03 
0.02 
0.01 
0.08 
0.05 


0.02 


0.01 
0.01 
0.02 
0.01 
0.05 
0.03 


0.04 
0.03 
0.04 
0.12 
0.05 
0.06 
0.20 
0.27 


HeMS 


0.01 
0.16 
0.01 
0.03 
0.02 
0.15 
0.04 
0.02 


0.03 


0.03 


0.07 


0.07 


HeHG 


0.08 


0.08 
0.09 
0.09 
0.01 
0.04 
0.08 


0.01 


HeGB 


0.01 


HeWD 


0.32 
0.38 
0.34 
0.33 
0.22 
0.33 
0.17 
0.23 


0.21 
0.47 
0.21 
0.22 
0.19 
0.42 
0.13 
0.16 


0.25 
0.56 
0.27 
0.35 
0.14 
0.33 
0.07 
0.11 


COWD 


0.09 
0.10 
0.10 
0.09 
0.08 
0.10 
0.07 
0.08 


0.07 
0.15 
0.07 
0.07 
0.08 
0.15 
0.05 
0.06 


0.01 
0.03 
0.01 
0.02 
0.01 
0.03 
0.01 
0.01 


HybWD 


0.37 
0.32 
0.39 
0.35 
0.42 
0.32 
0.37 
0.30 


0.17 
0.29 
0.16 
0.17 
0.21 
0.30 
0.13 
0.15 


0.13 
0.25 
0.14 
0.18 
0.18 
0.31 
0.10 
0.12 


Note. Companions (donors) in ULXs formed in models with constant SFR. Fractions for different companion types are provided separately for BHULXs and NSULXs. The age of the population is 10 Gyr. 
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Companions (Burst SFR 100 Myr Ago) 


Table 12 


Model 


ADO BN 
ADO BO1 
ADO BK 
ADO BS 
AD1 BN 
AD1 BOL 
AD1 BK 
AD1 BS 


ADO BN 
ADO BO1 
ADO BK 
ADO BS 
AD1 BN 
AD1 BOL 
AD1 BK 
AD1 BS 


ADO BN 
ADO BO1 
ADO BK 
ADO BS 
AD1 BN 
AD1 BOL 
AD1 BK 
AD1 BS 


MS 


0.89 
0.96 
0.97 
0.97 
0.90 
0.97 
0.98 
0.97 


0.86 
0.97 
0.94 
0.95 
0.85 
0.97 
0.94 
0.95 


0.81 
0.91 
0.86 
0.90 
0.79 
0.92 
0.87 
0.92 


HG 


0.11 
0.03 
0.03 
0.03 
0.10 
0.03 
0.02 
0.03 


0.11 
0.03 
0.04 
0.03 
0.12 
0.03 
0.03 
0.04 


0.05 
0.02 
0.01 
0.02 
0.06 
0.02 
0.02 
0.02 


RG 


CHeB 


0.02 


0.01 
0.01 
0.03 
0.01 
0.02 
0.01 


0.13 
0.06 
0.12 
0.08 
0.14 
0.05 
0.11 
0.06 


HeMS 


0.01 


0.01 


0.01 


0.01 


HeWD 


HybWD 


MS 


0.42 
0.01 
0.29 
0.07 
0.28 
0.14 
0.89 
0.73 


0.98 
0.97 
0.98 
0.96 
0.94 
0.99 
1.00 
0.99 


0.97 
0.88 
0.98 
0.85 
0.60 
0.98 
1.00 
0.98 


HG 
Zo 


Zo/10 


Z.,/100 


RG 


HeMS 


0.02 
0.98 
0.04 
0.67 


0.85 


0.10 
0.01 
0.12 


0.01 


HeHG 


0.49 
0.01 
0.57 
0.22 
0.59 
0.01 
0.10 
0.15 


0.01 
0.01 
0.01 
0.01 
0.04 


0.02 
0.01 
0.01 
0.01 
0.27 
0.01 


0.02 


HeGB 


0.01 


0.12 


0.01 


0.02 


0.13 


HeWD 


COWD 


HybWD 


0.06 


0.10 
0.03 


0.01 
0.01 


Note. Companions (donors) in ULXs formed in models with burst SF. Start of the burst was 100 Myr ago and its duration was 100 Myr. 
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Table 13 
Companions (Burst SFR 1 Gyr Ago) 


Model 


ADO BN 
ADO BO1 
ADO BK 
ADO BS 
AD1 BN 
AD1 BOL 
AD1 BK 
AD1 BS 


ADO BN 
ADO BO1 
ADO BK 
ADO BS 
AD1 BN 
AD1 BOL 
AD1 BK 
AD1 BS 


ADO BN 
ADO BO1 
ADO BK 
ADO BS 
AD1 BN 
AD1 BOL 
AD1 BK 
AD1 BS 


MS 


0.28 


0.51 


0.40 
0.32 
0.41 
0.40 
0.60 
0.64 
0.60 
0.58 


0.12 


HG 


1.00 
0.43 
1.00 
1.00 


0.01 


0.37 
0.44 
0.34 
0.37 
0.10 
0.22 
0.10 
0.12 


0.87 
0.47 
1.00 
0.90 


RG 


0.29 


1.00 
0.47 
1.00 
1.00 


0.18 
0.23 
0.19 
0.18 
0.24 
0.13 
0.24 
0.25 


0.03 
0.38 


0.02 


CHeB HeMS HewD HybWD 
0.01 
0.10 
0.03 
0.07 


MS 


0.02 
0.07 
0.01 
0.02 
0.01 
0.15 
0.03 
0.04 


0.01 
0.03 


0.01 
0.04 
0.08 
0.03 
0.02 


0.01 


0.19 
0.12 
0.05 


HG 
Lo 


(0) 
0.07 
0.03 
0.07 
0.12 
0.31 
0.17 
0.41 
0.47 

Zo/10 
0.74 
0.08 
0.78 
0.58 
0.47 
0.11 
0.73 
0.63 

Zo /100 
0.79 
0.10 
0.79 
0.60 
0.46 
0.08 
0.57 
0.46 


RG 


0.02 


0.01 
0.01 


0.01 


0.02 
0.02 
0.10 
0.12 


HeMS 


0.01 


0.01 


HeHG 


HeGB 


HeWD 


0.49 
0.12 
0.49 
0.46 
0.35 
0.09 
0.29 
0.26 


0.17 
0.32 
0.15 
0.28 
0.33 
0.28 
0.15 
0.23 


0.17 
0.36 
0.17 
0.32 
0.43 
0.29 
0.17 
0.30 


COWD 


0.09 
0.03 
0.09 
0.08 
0.07 
0.02 
0.06 
0.05 


0.03 
0.05 
0.03 
0.05 
0.06 
0.05 
0.03 
0.04 


0.01 
0.01 
0.01 
0.02 
0.03 
0.01 
0.01 
0.02 


HybWD 


0.33 
0.74 
0.33 
0.31 
0.25 
0.54 
0.20 
0.17 


0.05 
0.51 
0.05 
0.08 
0.11 
0.47 
0.05 
0.07 


0.03 
0.51 
0.03 
0.05 
0.07 
0.42 
0.03 
0.04 


Note. Companions (donors) in ULXs formed in models with burst SF. Start of the burst was 1 Gyr ago and its duration was 100 Myr. 
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Table 14 
Companions (Burst SFR 5 Gyr Ago) 
Model = MS HG RG CHeB. HeMS HeWD  HybWD MS HG RG HeMS  HeHG  HeGB  HeWD  COWD _ HybWD 
Zo 
ADO BN 0.37 0.63 
ADO BOl 0.01 0.75 0.24 
ADO BK 0.32 0.68 
ADO BS 0.19 0.81 
AD1 BN 0.31 0.69 
AD1 Bol 0.02 0.63 0.35 
AD1 BK 0.78 0.22 
AD1 BS 0.61 0.39 
Z5/10 
ADO BN 0.16 0.56 0.28 
ADO BOl 0.01 0.02 0.70 0.27 
ADO BK 0.17 0.55 0.28 
ADO BS 0.21 0.53 0.26 
AD1 BN 0.02 0.43 0.54 
AD1 Bol 0.01 0.04 0.54 0.41 
AD1 BK 0.25 0.34 0.41 
AD1 BS 0.38 0.28 0.34 
Zo /100 
ADO BN 0.10 0.25 0.66 
ADO BOl1 1.00 0.02 0.03 0.17 0.77 
ADO BK 0.04 0.30 0.66 
ADO BS 0.02 0.60 0.38 
AD1 BN 1.00 
AD1 Bol 0.17 0.09 0.10 0.68 0.08 
AD1 BK 0.38 0.62 
AD1 BS 0.03 0.97 


Note. Companions (donors) in ULXs formed in models with burst SF. Start of the burst was 5 Gyr ago and its duration was 100 Myr. 
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Model 


ADO BN 
ADO BO1 
ADO BK 
ADO BS 
AD1 BN 
AD1 BOL 
AD1 BK 
AD1 BS 


ADO BN 
ADO BO1 
ADO BK 
ADO BS 
AD1 BN 
AD1 BOL 
AD1 BK 
AD1 BS 


ADO BN 
ADO BO1 
ADO BK 
ADO BS 
AD1 BN 
AD1 BOL 
AD1 BK 
AD1 BS 


Note. Companions (donors) in ULXs formed in models with burst SF. Start of the burst was 10 Gyr ago and its duration was 100 Myr. A total absence of BHULXs may be observed. 


MS 


HG 


RG 


CHeB 


HeMS 


Table 15 
Companions (Burst SFR 10 Gyr Ago) 

HeWD HybWD MS HG 
Zo 

Z./10 

Zo/100 
0.05 0.07 
0.11 0.13 


RG 


0.40 
0.08 
0.38 
0.40 
0.19 
0.31 
0.92 
0.91 


0.07 


0.12 


1.00 
0.27 
1.00 
1.00 


HeMS 


HeHG 


HeGB 


HeWD 


0.60 
0.92 
0.62 
0.60 
0.81 
0.69 
0.08 
0.09 


0.83 
0.73 
0.83 
0.83 
0.89 
0.56 
0.90 
0.89 


1.00 
0.54 
1.00 
1.00 


0.01 


COWD 


0.17 
0.20 
0.17 
0.17 
0.11 
0.34 
0.10 
0.11 


HybWD 
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